Coherent generation of nonclassical light on a chip via photon-induced tunnehng and 

blockade 
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We report the observation of nonclassical light generated via photon blockade in a photonic crystal 
cavity with a strongly coupled quantum dot. By tuning the frequency of the probe laser with respect 
to the cavity and quantum dot resonance we can probe the system in either photon blockade or 
photon-induced tunneling regime. The transition from one regime to the other is confirmed by the 
measurement of the second order correlation that changes from anti-bunching to bunching. 

PACS numbers: 



Quantum dots in photonic crystals are interesting for 
their potential in quantum information processing [H |2] 
and as a testbed for cavity quantum electrodynam- 
ics. Recent advances in controlling [3^ and coherent 
probing m |5] of such systems open the possibility of re- 
alizing quantum networks originally proposed for atomic 
systems[6, 7, 8 . We recently showed[4j that the intensity 
of a probe beam coupled to a photonic crystal cavity is 
controlled by the presence of a single strongly coupled 
quantum dot that splits the cavity resonance into two 
eigenstates[9l llOj. Here we analyze the photon statistics 
of this probe beam. We show that the capture of a single 
photon into the cavity affects the probability that a sec- 
ond photon is admitted. This probability drops when the 
probe is positioned at one of the two energy eigenstates, 
resulting in photon antibunching. This result is analo- 
gous to the recent report on photon blockade of a neutral 
atom strongly coupled to a Fabry-Perot cavity [TT][T2]. In 
addition, we show that when the probe is positioned in 
between the energy eigenstates the probability of admit- 
ting subsequent photons increases, resulting in photon 
bunching. We call this process photon-induced tunnel- 
ing. This system represents an ultimate limit for solid- 
state nonlinear optics at the single photon level. Along 
with demonstrating the generation of nonclassical photon 
states, we propose an implementation of a single photon 
transistor [13] in this system. 

The optical system consists of a self-assembled InAs 
quantum dot (QD) with decay rate 7/27r — 0.1 GHz, cou- 
pled to a three-hole defect cavity [H] in a two-dimensional 
GaAs photonic crystal, as described in ref. |4]. The cou- 
pling rate g/lir = 16 GHz equals the cavity field decay 
rate K/27r = 16 GHz (corresponding to a cavity qual- 
ity factor Q=10,000), putting the system in the strong 
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coupling regime [HI HH]. We first characterize the sys- 
tem in photoluminescence (PL) by pumping the struc- 
ture above the GaAs bandgap. The PL scans in Figjl|b) 
show the anticrossing characteristic of strong coupling 
between the QD and the cavity. Here, the QD is tuned 
into resonance using local temperature tuning [15' around 
an average temperature of 20 K maintained in a contin- 
uous flow cryostat. We probe the system with linearly 
polarized laser beams (FigjlVa)) and observe the cross- 
polarized output as described in our previous work['I|. 
The measurement on the reflected port from this single- 
sided cavity is analogous to a transmission measurement 
in a Fabry Perot arrangement. 

The energy eigenstates of a two-level system strongly 
coupled to an optical resonator are grouped into two- 
level manifolds denoted |±,n), with energies hujn.± = 
h(niiJo±gy/n) , where n is the number of energy quanta in 
the system and wq is the bare cavity frequency (Figpla)). 
The anharmonic energy level spacing causes phenomena 
such as photon blockade |llj or photon-induced tunnel- 
ing. To observe photon blockade, a coherent probe beam 
(frequency Up) tuned at uji^± — loq ± g is coupled to the 
cavity. This probe is resonant with the first-order man- 
ifold, but detuned from transitions to the second mani- 
fold, (jJi^2 = tjQ ± g{V^ T 1) as shown in FigMa). Con- 
sequently, once a photon is coupled into the system, it 
suppresses the probability of coupling a second photon 
with the same frequency. As a result, the output field 
acquires sub-poissonian statistics with second-order cor- 
relation g^'^'{0) < 1, as recently demonstrated by Birn- 
baumm et al|llj in an experiment with neutral atoms. 

In addition to photon blockade, photon-induced tun- 
neling is expected near the bare cavity resonance {ujp — 
ujq — Aujp — !■ 0): the absorption of a first photon 
enhances the absorption of subsequent photons so the 
output consists of "photon bunches." In FigJ2'b) we 
show the output spectrum as the probe is tuned through 
the cavity and mark the resonance of the transitions 




FIG. 1: (a) Experimental setup. Laser pulses (40ps FWHM) are 
reflected from a photonic crystal cavity linearly polarized at 45° to 
the input polarization set by the polarizing beam splitter (PBS). 
The output light, observed in cross-polarization and carrying the 
cavity-coupled signal, was analyzed using a Hanbury-Brown-Twiss 
(HBT) setup that measures second-order correlation. The inset 
shows the suspended structure with the photonic crystal cavity and 
the metal pad for local temperature tunine[T3). (b) Anticrossing 
observed in photoluminescence as the QD is tuned into resonance 
with the cavity. The temperature tuning is done by linearly in- 
creasing the power P of the heating laser llSl . The red lines mark 
the cavity and QD resonance if they were decoupled. The right 
panel shows the spectrum at the anticrossing point marked by the 
blue line. 



|0) -^ \1,+) and |1,+) ^ |2,+). The simulated driv- 
ing field injects an average photon number (n) ~ 0.4 
when resonant with the polaritons in the first manifold, 
thus causing a slight saturation of the QD dipole. As 
AcOp —f 0, the probability of absorbing the first photon 
decreases. However, if a photon is nevertheless absorbed, 
it enhances the probability of capturing the second pho- 
ton. This process produces a photon-bunched output. 
This intuitive explanation can be extended in the same 
way to include higher order manifolds. 

The expected second-order correlation function for our 
system is shown in FigMc), where we plot the depen- 
dence of g'^''(0) for different detunings Aujp of the probe 
from the anticrossing frequency ujq . As expected from the 
intuitive argument above, the simulation predicts pho- 
ton bunching as Aup -^ 0. Photon blockade is evident 
in the antibunched region near Aup ^ il.Sg (see inset 
of FigMc)). However, the blockade does not occur for 
AtOp — ±g as previously explained, but at Aup ~ ±1.5g. 
The reason is that the linewidth of the eigenstates (^ k) 
is comparable to the splitting of the manifolds (~ 2g), 
resulting in significant overlap of the allowed transitions 



FIG. 2: (a) Energy diagram showing the flrst and second-order 
manifolds of the strongly coupled cavity/QD system. The en- 
ergy difference between consecutive manifolds is not constant as 
shown by the blue and the red arrows. This anharmonic spac- 
ing of the levels causes phenomena such as photon blockade llll 
and photon induced tunneling, (b) Simulated output intensity for 
a probe frequency tuned through the cavity/QD strongly coupled 
system. The red and blue lines indicate the transitions |0} — > |1, +) 
and |1,-|-} — > |2,+}. Photon blockade is expected at Aujp/g ^ 1 
because the absorption of a photon into |1, -I-) suppresses the prob- 
ability of absorbing a second photon of the same energy for a tran- 
sition to |2,+}. As Auip — > the absorption of a photon into the 
first manifold enhances the absorption probability into higher or- 
der manifolds, (c) Computed second-order correlation </"' (0) for a 
coherent laser probe reflected from the cavity. Inset shows that pho- 
ton blockade is expected when the probe detuning is Aup/g ~ 1.5 
(and not Aujp/g ~ 1) because of the finite linewidth of the po- 
laritons. As Acjp — > the output field is bunched, (d) Simulated 



time dependence of the second-order correlation for Aujp = 0. 
value for (g"' (r)) drops rapidly for time delays greater than ^ 
corresponding to the cavity photon lifetime. 
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between consecutive manifolds. 

We measure the time-dependent autocorrelation 
g^'^'ir) using the Hanbury-Brown-Twiss (HBT) setup 
shown in Fig [IT a) and described in ref . [ini HTj . The rele- 
vant features occur at the QD-cavity coupling rate g, en- 
veloped by the coherence time|18l. as shown in FigJ2jd). 
The coherence time for our system is given by the cav- 
ity photon lifetime 1/2k. Hence the time-dependent 
features in g^'^'ir) occur much faster than the 300 ps 
time-resolution of the single photon counting modules 
in the HBT setup. For this reason, we sample the au- 
tocorrelation function by short pulses {AtpwHM "^ 40 
ps, Alupwhm/'^t^ ~ 12 GHz) with a repetition rate of 
12.5ns. This probe pulse duration forms a compromise 
between fast sampling and a spectral linewidth narrow 
enough to resolve the relevant spectral features. In the 
remainder of the paper we present the measurements of 
g^^^r) for different detunings of the probe beam, denoted 
as g^^^T^Aujp/g). 

To observe photon blockade and photon-induced tun- 
neling, we measured the second-order correlation at de- 
tunings AcUp/g = —1.5 and Ausp/g = as shown Figls] 
The expected photon antibunching and bunching behav- 
ior is clearly visible at zero time delay (Fig|3lb,d)). The 
histograms also show bunching over timescales of hun- 
dreds of nanoseconds. This bunching results not from 
the quantum nature of the system, but from blinking be- 
havior of the QD. As reported by Santori et al[19j, such 
blinking behavior between a bright and a dark state re- 
sults in bunching near t=0, with a characteristic fall-off 
rate given by the mean switching rate. Our observations 
indicate that the blinking rates vary for different QDs and 
for different probe powers, and the quantum dot spends 
~ 80% of the time in the bright state. 

Photon blockade and photon-induced tunneling are 
quantified by the normalized second-order correlation 
function ^'•■^'(t, Awp/g). Each peak in the histogram of 
Figlslrepresents the unnormalized value of g^^^ (r, Aujp/g) 
averaged over the pulse duration of 40 ps. We express 
this time averaging by using the notation ^"'(t, AcUp/g). 
Because of the QD blinking there are two choices for the 
normalization constant. One choice is to normalize such 
that 5*-^-'(t -^ oo,Aujp/g) = 1. We keep the notation 
^"'(t, Awp/g) for this normalized quantity (FiglSTb)). 
We stress that g^^> (r, Aujp/g) captures both the quan- 
tum and classical nature (i.e. blinking) of the output 



field. The normalization constant Gqc for g^^'{T,Aujp/g) 
is determined by fitting the histogram with the func- 

tion G^^\mTo) = (G^'^ - G^^xp [-mTo/T] + G^^ for 

(2) 

m > 1. The quantity G (t) represents the number of 
counts at time ttiTq, where m indexes the peak number 
with TO = corresponding to t = 0, and Tq = 12.5 ns 
is the pulse repetition period. The other choice of nor- 

(2) 

malization constant is Gq , which is equivalent to nor- 
malizing to the nearest neighbor peaks at t = ±12.5 ns. 
When referring to this normalization we use the nota- 



tion (7q '{T,Aujp/g) (Figlstb)). The time-dependence in 
g'^'^\T^Au}p/g) on the time scale of 12.5 ns is due almost 
exclusively to the fast interaction between probe and QD 
and is only negligibly affected by the comparatively slow 
QD blinking process. Thus the nearest-neighbor nor- 
malized ^q (r, Aojp/g) better captures the photon block- 
ade (^Q (0, Acjp/g) < 1) and photon-induced tunneling 
(g('^(0,Aa;p/5)>l) effects. 

In the case of photon blockade, g^^^O, 1.5) = 0.912 ± 

0.005 and g^^\Q,l.b) = 0.881 ± 0.009 (Figjajb)), show- 
ing the antibunched quantum nature of the system for 
both choices of normalization. For photon-induced tun- 
neling, g(2)(o, 0) = 1.33±0.02 and 51,^^(0,0) = 1.27±0.02 
(FigMd)). These values are different from the ideal 
theoretical prediction in FigJ2Fc) because of the back- 
ground noise caused by the imperfect extinction of the 
cross-polarized experimental setup (signal to noise ratio 
^6:1), QD blinking and finite bandwidth of the probe. 
Both the background and the output signal when the QD 
is in the dark state have g*^^\T,Au)p/g) — 1. 

We repeated the autocorrelation measurements for 
a large set of detunings to map the full spectrum of 

(2) 

5q (0, Awp/g). The measurement of the full autocor- 
relation spectrum entails several challenges such as sam- 
ple drift resulting in fluctuating coupling intensity into 
the cavity, and fluctuating temperature. To map the de- 

(2) 

pendence of TJq {0,Aujp/g) on probe detuning, we took 
care to maintain constant coupling into the cavity mode 
for the full duration of the experiment. Our most con- 
sistent data set is presented in Fig|4] where we plot 

(2) ^"^ 

5g {0,Aojp/g) for different detunings of the probe fre- 
quency. We kept a constant probe power of ~ 1.0 nW 
before the objective lens corresponding to an average cav- 
ity photon number (n) ^ 0.4, and the coupling was re- 
optimized for every data point. The lowest value for 

(2) 

5q {0,Aujp/g) obtained in this data set is not as an- 

tibunched as the value g^ (0,1.5) — 0.88 reported in 
FigJ3[b), mainly because we could not reproduce exactly 
the same coupling conditions. We found that the data 
in Fig|4] is well fitted by a numerical model that takes 
into account pulses of finite bandwidth, QD blinking and 
background from the imperfect extinction of the cross- 
polarized setup (see Methods section for details). 

The experimental data in Fig|4] shows that this 
strongly coupled system allows control of the statistics of 
the output field from sub-poissonian to super-poissonian. 
Thus, the probability of transmitting a specific number 
of photons can be tailored by controlling Aujp. Given 
the above-mentioned properties, this device can be used 
as a single photon transistor [ISj where the presence of a 
gate field could control the transmission of a signal field. 
One scheme to implement the transistor is to tune the 
frequency of the gate field resonant with one of the po- 
laritons in the first-order manifold, say u;o + g. A photon 
injected at loq + g, increases the probability of absorbing 
photons that are resonant with the |l,-l-) -^ |2, -I-) tran- 
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FIG. 3: Measurement of the second-order correlation function for coherent laser pulses reflected from the photonic crystal cavity with 
a strongly coupled QD. (a) Photon blockade, manifested here in the sub-poissonian statistics, observed when the probe is detuned by 
Auip/g = 1.5. (b) Beside the antibunching at r = the data in panel a shows bunching due to QD blinking over timescales of hundreds 
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7^(2) 



of nanoseconds. The red and blue lines show the normalization levels G^o and Gg observed as t — > oo and as r — > respectively. 
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The second-order correlation is g'^^'{Q, 1.5) = 0.91 when normalized to G)^' and Qq'' {0, 1.5) = 0.88 when normalized to Gq'' (c) Photon- 
induced tunneling is observed when the laser pulse is tuned at the anticrossing wavelength Aujp/g = 0. (d) Depending on the choice of 

normalization for the data in panel c, the second-order correlation takes the values g"'(0,0) = 1.34 and Qq (0,0) = 1.27. 




FIG. 4: Second-order correlation function g^ ' {0, Aujp/g) for dif- 
ferent detunings between the probe and the anticrossing frequency. 
As the probe is tuned through the resonance of the QD-cavity sys- 
tem, the output field shows antibunched and bunched behavior as 
it transitions from photon blockade regime to photon-induced tun- 
neling regime. The fit takes into account the finite pulse-width of 
the probe, QD blinking and background due to the imperfect ex- 
tinction of the cross-polarized setup. In the bottom-right corner 
we show the relative width of the pulsed laser probe. 



sition at wq + g{V^ ~ !)• If the signal is tuned to this 
frequency, the presence of the gate field enhances the 
transmission of the signal field. The photonic crystal ar- 
chitecture allows for easy integration of such a single pho- 
ton transistor with photonic crystal waveguides [^Hl HI] so 
the single photon switching is done directly on the chip. 
The most straightforward configuration would be a pho- 
tonic crystal cavity butt-coupled in between two photonic 
crystal waveguides j22j. 

In conclusion we observed the phenomena of photon 
blockade and photon-induced tunneling in a solid-state 
strongly coupled QD/cavity system. By changing the 
frequency of the input beam, we synthesized nonclassical 
states of light with various degrees of bunching or anti- 
bunching (0.91 < g'^^O) < 1-34). These results comple- 
ment the photon blockade experiments in atomic physics 
by showing the full dependence of g^'^'ir) on the probe 
frequency. They have promise in realizing single pho- 
ton transistors for on-chip optical logic devices operat- 
ing at the single photon level, or the implementation of 
two-qubit gates and generation of nonclassical states for 
quantum information processing, quantum lithography 
and quantum metrology. 



Methods 



Autocorrelation Measurement: We scan several cavi- 



ties until we find one which contains a strongly coupled QD, 
as determined by the anticrossing behavior in photolumines- 
cence between QD and cavity during temperature tuning. 
Then we direct the pulsed laser beam at the cavity and ob- 
serve the reflected beam in cross-polarization. While tuning 
the local temperature with an additional heating beam, we 
adjust the probe beam coupling to optimize the QD-induced 
reflectivity drop, as described for the continuous-wave beam 
in ref. [4]. Then we stop scanning and temperature-tune 
the QD and cavity onto resonance. With the pulsed probe 
beam at different detunings with respect to the anticrossing 
point, we measure the autocorrelation signal by passing the 
reflected probe through a grating filter (to remove stray 
light) followed by the HBT setup. To limit sample drift, the 
alignment procedure is repeated for every data point in Fig. 4. 

Data Analysis: The numerical model for the second-order 
coherence in Fig. 4 is based on numerical integration of the 
quantum master equation. A time-dependent driving term in 
the Hamiltonian represents the 40-ps excitation pulses. The 
intensity of the drive field matches the intensity used in the 
experiment, representing one-third of saturation. In our ex- 
periment, this intensity was 1 nW for the incident beam, mea- 
sured before the objective lens. The state of the QD/cavity is 
time-evolved using a quantum Monte Carlo approach, which 
we based on the qotoolbox of ref. [23]. The QD can emit into 
free space or into the cavity mode, which in turn dissipates 
energy into the output channel at the loss rate co/Q. We 
then compute the autocorrelation on the output channel. The 
simulation additionally accounts for QD blinking and laser 
background. The full second-order coherence is calculated as 
weighted sum of the different contributions. 



G^^Ur) 
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(1) 



PBG'^>(T)+pBGG'-^>ai^)+PDG'^>{T) (2) 



^(2), 



(2), 



where the autocorrelation function Gg (t) accounts for the 
QD bright state, G]-, (r) for the QD-dark state (calculated 
using g — > 0), and GgQ{T) for background laser signal (a 
coherent state). We then normalize by the simulated nearest 

neighbor peaks (Go ) to obtain Qq (t), as shown for example 
in Fig[3j;d). 
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